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Adsorpcja na powierzchniach heterogenicznych. Metody ilościowej oceny 
rozkładu energii adsorpcji na powierzchni adsorbentów rzeczywistych 
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Adsorption on Heterogeneous Surfaces. Quantitative Estimation of the Distribution 
of Adsorption Energy on Real Surfaces 


Do ilościowej charakterystyki energetycznej heterogeniczności po- 
wierzchni adsorbentu przyjęto w teorii adsorpcji tzw. różniczkowy roz- 
kład energii adsorpcji. Rozkład ten oznaczany zwykle przez y(e) definiuje 


się jako A , gdzie w jest ilością centrów adsorpcyjnych na powierzchni 


adsorbentu. Z matematycznego punktu widzenia wyznaczanie funkcji roz- 
kładu y(e) sprowadza sie do rozwiązania równania calkowego, definiują- 
cego izotermę globalną ©;+(p) na płatowo heterogenicznych powierzch- 
niach [1], 


PIN Nn [à lp E) X: [E) d W 


W powyższym równaniu 6;(p, e) jest izotermą lokalną, wyrażoną w jed- 
nostkach względnego pokrycia, N+(p) jest ilością adsorbatu w fazie po- 
wierzchniowej, Nm jest pojemnością monowarstwy, natomiast ( jest ob- 
szarem możliwych zmian energii adsorpcji e. Dalej przez y(sc) będziemy ro- 
zumieć funkcję normalizowalną do jedności: 


Ją y(s)de — 1. 
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W poprzedniej pracy [1] dyskutowano analityczne rozwiązanie równa- 
nia (1) ze względu na izotermę globalną ©;(p) lub funkcję rozkładu energii 
x(e), przy założonych a priori analitycznych wyrażeniach dla dwu pozosta- 
łych funkcji występujących w równaniu (1). Pierwszymi fundamentalny- 
mi pracami na ten temat są prace Sipsa [2, 3], sprowadzajace problem 
wyznaczania funkcji (e) do wyliczenia odwrotnej transformaty Stieltjesa. 
Istotny postęp w zakresie wyliczania funkcji y(e) z danych eksperymen- 
talnych nastąpił dopiero w kilkanaście lat po ukazaniu się prac Sipsa, kie- 
dy to pojawiło się szereg nowych, numerycznych metod wyznaczania roz- 
kładu energii adsorpcji. 

Celem tej pracy jest krytyczna dyskusja numerycznych metod wyzna- 
czania funkcji x(e) z eksperymentalnych izoterm adsorpcji. Jest to problem 
obecnie szczególnie interesujący, biorąc pod uwagę szerokie wykorzysta- 
nie elektronicznych maszyn cyfrowych we współczesnych badaniach nau- 
kowych. Dyskutowane poniżej metody numeryczne mają charakter ogólny 
i mogą być użyteczne do wyznaczania wielu wielkości fizykochemicznych, 
takich jak: nadmiar powierzchniowy w adsorpcji z roztworów [4], izoter- 
my adsorpcji z mieszanin gazowych [5, 6], kinetyczna izoterma adsorpcji 
na heterogenicznych powierzchniach [7], drugi współczynnik wirialny dla 
heterogenicznych powierzchni [7], pole powierzchni i porowatość ciał sta- 
łych [8], globalna objętość retencji [9], emisja elektronów z powierzchni 
metali [10], funkcja rozkładu energii w biomolekułach [11], funkcja roz- 
kładu czasu relaksacji charakteryzująca zachowanie elektryczne dielek- 
tryków [12] oraz rozkład odległości donor-akceptor [13]. Ponadto wspom- 
niane metody okazują się pomocne w dyskusji efektów cieplnych promie- 
niowania [14], reakcji antygen—antyciało [15], radiograficznej lokalizacji 
raka [16] oraz w radiometrycznych pomiarach składu atmosfery [17]. Do 
opisu wymienionych zjawisk stosuje się wyrażenia analogiczne do równa- 
nia (1). W dalszej części pracy ograniczymy się do przedstawienia naj- 
ważniejszych metod wyznaczania funkcji rozkładu energii adsorpcji y(e), 
znajdujących szerokie zastosowanie do ilościowej charakterystyki ukła- 
dów adsorpcyjnych. 


METODA ADAMSONA 


Metoda Adamsona [18] sukcesywnych graficznych lub ogólniej — 
numerycznych aproksymacji jest dalszym rozwinięciem idei przybliżenia 
kondensacyjnego, wprowadzonej przez Rogińskiego [19]. Rozważa- 
nia prowadzone są tu w terminach całkowego rozkładu energii X(e) okreś- 


lonego równaniem: 
X(E) KIE | Q) 
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Podstawowe calkowe równanie (1) dla globalnego efektu adsorpcji na po- 
wierzchni heterogenicznej można zapisać w następującej, równoważnej 
formie: 


4 
8,ip)= [8,(p,E)dX (3) 
0 


Jako pierwsze przybliżenie dla funkcji X(e) przyjmuje sie funkcję otrzy- 
maną w metodzie przybliżenia kondensacyjnego. Procedura otrzymywania 
dalszych przybliżeń jest zbliżona do metody Berenyi [20] określania 
krzywych charakterystycznych w teorii potencjału. 

Przyjęcie określonej izotermy lokalnej pozwala przedstawić ©, jako 
funkcję e dla każdej zadanej wartości ciśnienia. Pierwsze przybliżenie X 
znajduje się z eksperymentalnej izotermy adsorpcji O;. Dalej, wykorzy- 
stując krzywe X0— i O©r—e, sporządza się dodatkową krzywą zależności 
©, względem odpowiadającej jej wartości XU. Pole powierzchni pod tą 
krzywą daje obliczoną wartość ©, (oznaczmy: 0,7%) odpowiadającą zada- 
nej wartości ciśnienia p. Postępowanie to powtarza się dla wszystkich 
przyjętych wartości ciśnienia. Wartości 0,7%! są w najogólniejszym przy- 
padku dość różne od wartości eksperymentalnych 0,***5. Dla każdej war- 
tości p drugie przybliżenie funkcji X(e) oblicza się z następującej zależ- 
ności: 


eks 


x 6) = X"(E): d (4) 


W celu znalezienia trzeciego przybliżenia X? powtarza sie powyższą pro- 
cedurę, używając wykresu funkcji X© względem energii e. 

Szczegóły metody Adamsona zmieniają się znacznie, gdy zmienia 
się równanie opisujące adsorpcję lokalną. Użycie każdej nowej izotermy 
lokalnej jest w tej metodzie zabiegiem przeważnie trudnym. 

Najprostszym przykładem ilustrującym metodę Adamsona jest 
prawdopodobnie przypadek lokalnej izotermy Langmuira: 


6, (a£) = ap (1 * ag) 6) 
W metodzie tej wygodniej jest operować funkcjami z i 65, zależnymi od 


stałej a, która, jak wiadomo, wiąże się z energią adsorpcji e w następujący 
sposób: 


a= [aK exo [m (8) 
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W równaniu (6) K jest stala temperaturowo zalezna, zwiazana z moleku- 
larną funkcją rozdziału molekuł warstwy pierwszej [1]. 

To „langmuirowskie” wydanie metody Adamsona sprawdzono, 
obliczając y(e) dla różnych układów adsorpcyjnych. Do obliczeń wyko- 
rzystano między innymi dane Joynera i Emmetta [21] dla układu 
azot—Graphon w temp. 78?K. Z otrzymanego rozkładu energii wynika, że 
powierzchnia tego adsorbentu jest prawie homogeniczna. Porównanie cał- 
kowego rozkładu energii X(e) z rozkładem otrzymanym z pomiarów kalo- 
rymetrycznych wykazuje zadowalającą zgodność obu tych rezultatów. 

W dalszych badaniach Adamson i współprac. zmodyfikowali me- 
todę sukcesywnych aproksymacji tak, by można było w niej użyć równa- 
nie Fowlera-Guggenheima [22, 23] lub BET [24] w charakterze izotermy 
lokalnej. Aby sprawdzić poprawność swojej zmodyfikowanej metody w 
przypadku równania BET, Adamson i wspólprac. [22] wykonali obli- 
czenia funkcji rozkładu energii y(e) dla danych adsorpcyjnych Draina 
i Morrisona [25] opisujących adsorpcję argonu na rutylu w temp. 
85?K. Drain i Morrison, wykorzystując dane kalorymetryczne, 
określili funkcję x(e) dla tego układu. Przy pomocy przybliżenia konden- 
sacyjnego znaleźli funkcję x(e) dla szeregu temperatur, a następnie ekstra- 
polowali swoje wyniki do temperatury zera bezwzględnego. Porównanie 
wyników Adamsona i współprace. [22] z rozkładem energii znalezio- 
nym przez Draina i Morrisona [25] pokazuje, że obie krzywe roz- 
kładu energii są raczej podobne, chociaż funkcja Adamsona i współ- 
prac. jest przesunięta na skali energii (patrz ryc. 1). Według tych ostat- 


X(£jcm? moL/fcau:g 
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Ryc. 1. Funkcje rozkładu energii adsorpcji dla argonu na rutylu w temp. 85?K; 
1 — funkcja x*(e) wyznaczona przez Dormanta i Adamsona [24], 2 — funkcja 
Draina i Morrisona [25] 
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nich, najbardziej prawdopodobnego wytłumaczenia tego faktu należy do- 
szukiwać się w dowolnym założeniu teorii BET, że ciepło adsorpcji w dru- 
giej warstwie adsorpcyjnej jest równe ciepłu parowania. Innym możli- 
wym źródłem tego przesunięcia może być wybór BET-owskiej stałej c(e). 
Rozkłady wykazujące znacznie lepszą zgodność z rezultatami Draina 
i Morrisona można otrzymać przyjmując, że ciepło adsorpcji w drugiej 
warstwie jest o 150 cal/mol większe niż energia parowania, lub też zakła- 
dając, ze przedeksponencjalny faktor we wzorze na stałą c(e) jest nieznacz- 
nie większy od jedności. 

Metoda Adamsona jest pierwszą, szeroko akceptowaną, numerycz- 
ną metodą obliczania funkcji rozkładu energii adsorpcji z(e) z tabelarycz- 
nych 'danych eksperymentalnych [22, 24, 26]. W zasadzie metoda ta po- 
zwala na użycie dowolnej funkcji w charakterze izotermy lokalnej, cho- 
ciaż przejście od jednej izotermy lokalnej do innej nie jest procedurą 
prostą ze względu na konieczność znalezienia funkcji skokowej, aproksy- 
mującej jak najlepiej użytą izotermę lokalną. 


METODA VAN DONGENA 


Metoda ta polega na przybliżaniu funkcji rozkładu energii z(e) wy- 
kładniczym wielomianem postaci [27] 


x(€)=exp [> b; gt) (7) 


W wielomianie tym występuje m parametrów: by, ..., bm, umożliwiają- 
cych jak najlepsze dopasowanie go do rzeczywistej funkcji 4(e). Wybór 
przez Van Dongena właśnie takiej analitycznej aproksymacji był raczej 
intuicyjny. Okazuje się jednak, że ma on racjonalne, teoretyczne podsta- 
wy. Przyjmijmy wykładniczo-logarytmiczne rozwinięcie dla funkcji ©;(p), 
postulowane przez Heera [28], 


N, (p)= exp pa 8j (RT ln £y] (8) 
fz 


gdzie B; są współczynnikami rozwinięcia, a pm jest ciśnieniem charakte- 
rystycznym, Zgodnie z metodą przybliżenia kondensacyjnego takiej po- 
staci analitycznej dla globalnej izotermy adsorpcji ©, odpowiada nastę- 
pująca funkcja rozkładu [29, 30]: 


XE) [5 (6-6) ""Jeo(Sa & E] © 
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gdzie 


£,, * RT Un d oraz x (E) = 26 (E)N (10) 


Można zauważyć, że przybliżenie Van Dongena jest uproszczoną formą 
zupełnie ogólnego, analitycznego wyrażenia na funkcję (c). 

W numerycznej metodzie Van Dongena przyjmuje się, że granice cał- 
kowania w równaniu (1) są skończone. Jeśli p, i pg są odpowiednio dolną 
i górną granicą przedziału ciśnień, to e, i eg (granice całkowania w rów- 
naniu (1) określamy korzystając z metody przybliżenia kondensacyjne- 
go. W celu obliczenia parametrów b; Van Dongen posłużył się metodą naj- 
mniejszych kwadratów. Wyliczenia iteracyjne wymagają znajomości ciągu 
startowych wartości współczynników b;. Dla ich określenia Van Dongen 
zaproponował dwa różne, praktyczne sposoby postępowania. 

W pierwszym z nich punktem wyjściowym obliczeń jest stały rozkład 
energii adsorpcji, a więc b,(? jest równe In[N4,(c,—eg)], a wszystkie wyż- 
sze współczynniki są równe zeru. Druga procedura jest bardziej efektyw- 
na. Jako pierwsze przybliżenie funkcji y(e) przyjmuje się tu kondensa- 
cyjną funkcję rozkładu. Przybliżenie pierwszego rzędu dla współczynni- 
ków b, znajduje się przez numeryczną aproksymację krzywej kondensa- 
cyjnej wielomianem (7). Iteracyjna metoda Van Dongena daje szybko 
zbieżne rozwiązanie pod warunkiem, że: 

1) granice całkowania wybrane są zgodnie z badanym przedziałem 
ciśnień, 

2) stopień wielomianu aproksymującego jest nie większy niż 6, 

3) liczba punktów eksperymentalnych jest większa niż 30, 

4) punkty te odczytuje się z wyrównanej postaci izotermy. 

W charakterze izotermy lokalnej Van Dongen użył poliwarstwowej 
izotermy Broekhoffa [31], będącej uogólnieniem monowarstwowej 
izotermy Hilla-de Boera. W modelu adsorpcji, w oparciu o który wypro- 
wadzono równanie Broekhoffa, przyjmuje się, że faza zaadsorbowana skła- 
da się z szeregu oddzielnych warstw [32]. Molekuły zaadsorbowane w każ- 
dej z tych warstw zachowują się jak dwuwymiarowy gaz Van der Waalsa. 
Model ten uwzględnia wzajemne oddziaływania między molekułami ad- 
sorbatu oraz zanikanie potencjału adsorpcyjnego wraz ze wzrostem od- 
ległości od powierzchni. Jeśli chodzi o odległości między poszczególnymi 
warstwami, to przyjmuje się, że są one jednakowe i równają się odległości 
pierwszej warstwy adsorpcyjnej od powierzchni. Ostatnią odległość ocenia 
sie zwykle jako średnicę zderzenia molekuł adsorbatu. Van Dongen wpro- 
wadził do tego modelu jeszcze jedno uproszczenie zakładając, że częstość 
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wibracji w kierunku prostopadłym do powierzchni jest dla wszystkich 
warstw taka sama i nie zależy od stopnia pokrycia powierzchni. W ilustra- 
cyjnych obliczeniach numerycznych Van Dongen wykorzystał dane ekspe- 
rymentalne opisujące układy azot—aerosil i krypton—aerosil w temp. 
77,4°K, oraz azot—aeroal i krypton—aeroal w temp. 77,4?K. Funkcje roz- 
kładu energii otrzymane dla dwóch ostatnich układów zmierzają do nie- 
skończoności, gdy e dąży do zera. Wynik taki nie posiada oczywiście sensu 
fizycznego. Wydaje się, że w niektórych przypadkach sama funkcja (7) 
może być przyczyną otrzymywania takich kłopotliwych rozkładów. Apro- 
ksymacja rzeczywistej funkcji (e) wielomianem (7) nie jest efektywna, 
gdy z(e) składa się z wielu dobrze oddzielonych pików. Z tego punktu wi- 
dzenia metoda Van Dongena powinna być użyteczna, gdy y(e) posiada sze- 
reg blisko siebie leżących pików energetycznych. 

W tym miejscu należy wspomnieć o metodzie numerycznej Rossa 
i współprac. [33], dotyczącej charakterystyki adsorbentów wykazujących 
gaussowski rozkład energii, tj. równanie (7) dla b,7-2—0. Powyższa meto- 
da jest doskonale opracowana w monografii Rossa i Oliviera [33]. 


UOGÓLNIONA METODA TRANSFORMATY STIELTJESA 


Ostatnio opracowano dwie nowe, numeryczne metody obliczania (s) 
z danych eksperymentalnych, które można określić jako uogólnienie me- 
tody odwracania transformaty Stieltjesa. Autorami pierwszej z nich są 
Rudziński i współprac. [34, 35]. Rozważali oni najogólniejszy przy- 
padek modelu BET, gdy na powierzchni adsorbentu tworzy się skończona 
liczba (n) warstw adsorpcyjnych. Równanie izotermy ma w tym przypad- 
ku następującą postać: 


_awmy_,_ c(&):g 
8, (x€) 8 Unx 1* c(E):4 e 
gdzie: 


gu r 
J 

je az y=? X (12) 

X=P/ps2razj A 


W powyższych równaniach p, jest ciśnieniem pary nasyconej, natomiast 
n oznacza liczbę tworzących się warstw adsorpcyjnych. Zasadniczym eta- 
pem metody Rudzińskiego jest aproksymowanie eksperymentalnych da- 
nych adsorpcji następującym wyrażeniem: 


12 Annales, sectio AA, vol. XXXI/XXXII 
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ARK eiog T 1 FU 
0, (x) "dlny ^ 3 Lu (13) 


gdzie l jest pewną liczbą całkowitą. Obliczona z odwrócenia transformaty 
Stieltjesa funkcja x(e) dana jest wyrażeniem: 


4 ipergziące isi qq | Lg 
«-nhe- Lc eo [E oo [i] 
gdzie (14) 

od eor) 


Mozna zauważyć, że: 
X->0 (15) 


Zatem metoda Rudzińskiego polega na aproksymowaniu eksperymental- 
nej izotermy szeregiem izoterm Freundlicha. Do ilustracyjnych obliczeń 
numerycznych wykorzystano dane eksperymentalne dla układu argon—żel 
krzemionkowy w temp. 78?K otrzymane przez Aristowa i Kisie- 
lewa [36]. Występowanie dwóch maksimów na krzywej rozkładu energii 
związane jest prawdopodobnie z istnieniem różnych rodzajów grup hydro- 
ksylowych na powierzchni żelu krzemionkowego. 

Drugą metodą opracowaną na bazie transformaty Stieltjesa jest meto- 
da Jarońca [30]. W metodzie tej aproksymuje się izotermę adsorpcji 
globalnej 9, wyrażeniem (8). W obliczeniach praktycznych wykorzystuje 
się oczywiście pewną sumę częściową zamiast nieskończonego szeregu. Ja- 
roniec ogranicza swoje rozważania jedynie do adsorpcji monowarstwowej, 
przyjmując do opisu adsorpcji lokalnej równanie Langmuira. Sprowadza- 
jąc równanie całkowe (1) do transformaty Stieltjesa otrzymuje się funkcję 


rezkładu energii w następującej postaci: (16) 
4e. 1 A) d En- EW x (Ji 2 S! fen- p U-2 
xte): wm [ove 6). Enn) - (1). (Eng Y]. 
J” 
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W ilustracyjnych obliczeniach numerycznych Jaroniec wykorzystał dane 
eksperymentalne Aristowa i Kisielewa [36] dla układu azot—żel 
krzemionkowy w temp. 78?K. Funkcja rozkładu x*(e), charakteryzująca 
ten układ adsorpcyjny jest przedstawiona na ryc. 2. 


05 


2000 2500 3000 3500 
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Ryc. 2. Funkcje rozkładu energii adsorpcji dla azotu na żelu krzemionkowym 
w temp. 78?K; 1 — funkcja z*(e) wyliczona z równania (16), 2 — funkcja obliczona 
z równania (9) 


Porównanie dwóch przedstawionych wyżej numerycznych metod obli- 
czania xy(c) sprowadza się praktycznie do porównania sensu aproksymacji 
(13) i (8). Przybliżenie Jarońca można określić jako uogólnienie równania 
Dubinina-Raduszkiewicza [29]. Przy dostatecznie niskich 
wartościach ciśnień wszystkie układy adsorpcyjne powinny spełniać rów- 
nanie DR, niezależnie od tego, jaki jest ich dokładny rozkład energii. Tak 
więc aproksymacja (8) będzie nadzwyczaj efektywna, jeżeli będziemy dy- 
sponować danymi eksperymentalnymi z obszaru niskich pokryć [37]. Z 
kolei przy umiarkowanych i średnich ciśnieniach adsorbatu układy adsorp- 
cyjne powinny wykazywać tendencję do zachowywania się zgodnie z rów- 
naniem izotermy Freundlicha [37]. W tym przypadku aproksymo- 
wanie funkcji ©:+(p) szeregiem izoterm Freundlicha będzie skuteczniejsze. 

Użycie jednej z dwóch przedstawionych wyżej metod zależy więc od 
rodzaju danych eksperymentalnych, którymi dysponujemy do obliczania 
funkcji z(e). 
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METODA OPARTA NA ROZWIJANIU FUNKCJI PODCALKOWEJ 
W POJEDYNCZY I PODWOJNY SZEREG FOURIERA 


Ostatnio Rudziński i wspólprac. [38] zaproponowali nową, zupel- 
nie ogólną, numeryczną metodę rozwiązywania równania (1). Metoda ta 
opiera się na rozwijaniu funkcji ©,(p, e) w podwójny szereg Fouriera ze 
względu na p i e, a ż(e) w pojedynczy szereg Fouriera. Aby można ją było 
zastosować, muszą być spełnione następujące warunki: 


kim, & (p £)*1 oraz tim &(p) 1 | (18) 


Analityczne równania izoterm, opisujących adsorpcję monowarstwową, 
spełniają oczywiście warunek pierwszy. Warunek drugi jest najwyraźniej 
ograniczeniem, gdyż zwykle zakłada się, że najwyższym możliwym ciśnie- 
niem jest ciśnienie pary nasyconej ps. Założenie to przyjmuje się trady- 
cyjnie od wielu lat, chociaż model BET nie jest wcale ograniczony do 
ciśnień mniejszych niż ps. Ograniczenie to wynika z pewnego dodatkowe- 
go założenia, że faza zaadsorbowana zachowuje się jak faza ciekła przy 
p=ps. Dobrze. znane prace Andersona [39] oraz Brunauera, 
Skalnego i Bodora [40] wskazują na konsekwencje tego błędnego 
założenia. 

Równanie znanych multiwarstwowych izoterm O;V(p, e), np.: izoterm 
BET, Jovanovica, można zapisać w postaci iloczynu dwóch funkcji 6,(p, e) 
i p(p). Funkcja Oi(p, e) spełnia warunek pierwszy (18), natomiast funkcja 
p(p) opisuje efekt adsorpcji poliwarstwowej. W tej sytuacji iloraz multi- 
warstwowej adsorpcji © (p) przez funkcję g(p) spełnia warunek drugi 
(18). Zgodnie z powyższymi ustaleniami równanie całkowe (1) będzie miało 
postać: 


nlp) =N; (P/P (p) ef alp E) X (E)d6 (19) 


Funkcje N:(p) i ©,(p, e) czynią zadość warunkom (18), a więc warunki te 
nie stanowią żadnego ograniczenia dla tej metody. Rozwiązanie równania 
calkowego (1) otrzymuje sie przez rozwinięcie funkcji podeałkowych w 
szeregi Fouriera, stosując ortonormalny i zupełny układ funkcji u, okreś- 
lonych w przedziale [a, f]. Z teorii szeregów Fouriera wiadomo, że takie 
ortonormalne i zupełne układy funkcji istnieją w dowolnym przedziale 
[a, 8]. Dokonajmy następującego przekształcenia zmiennych: 


„€=fft) ; tefa] oraz p=9(5) ; sejal ] (20) 
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gdzie f(t) i g(s) są ciągłymi i jednoznacznymi funkcjami określonymi w 
przedziale [a, 8] tak, że 


fFla)-€ 6a F(P) =€g (21) 
9(a) =0 ; HP) *5, 
gdzie pm jest najwyższą wartością p, dla której znamy eksperymentalną 
wartość ©+(p). e. i eg oznaczają najmniejszą i największą wartość energii 
adsorpcji na danej powierzchni. Celem powyższych przekształceń jest za- 
miana przedziałów [e», eg] i [0, Pm] we wspólny przedział [a, 8]. Równa- 
nie całkowe (19) przyjmuje następującą formę: 


B ^ 
G()- fH (4 s) F[&) dt (22) 
a 
"^ — e" ^ 
H(t,s) = 6, (918) Ft) plt) | (24) 


X x, (F(&)) Plé) / pit) (25) 
Z fizycznego znaczenia Silinir F(t) i G(s) wynika, że: 


F(t) € - [ap] oraz Cí(s)e Leja, B] 


Nieznana funkcja q(t) określona jest tak, aby spełnione były warunki: 


BB 
ia Hilt, s) dt ds <= (26) 
ad, 
oraz 
88 2 
JJ [H (ć,s) ] déds < œ (27) 


Wiadomo, że jeśli f ver: „jest zupełnym układem ortonormalnych. funk- 
cji w przedziale [a, £], to iloczyn fui: (t) * we(s)| jest zupełnym i ortonor- 
malnym układem funkcji w przedziale [a, 8] X fa, 8]. Dlatego też uwzgled- 
niając warunek (27) H(t, s) można rozwinąć w podwójny szereg Fouriera 
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H (6,5) 23M u; (t): a, (s) . (28) 


L=0 k:0 
gdzie hi, są współczynnikami rozwinięcia danymi warunkiem: 
BB 
K. 29 
hik - f [^(&3)-; (€) u (s)d6ds i 
Aa A 
Podobnie 
oo 
G(s) = Du: DEDOS (30) 
k=0 : 
F(t) =), foru ltl (31) 
6*0 


Ponieważ funkcja G(s) jest znana, współczynniki gy oblicza sie z równania: 


e 
= A i 32 
| de] Gla) u (s)ds 62) 
o 
Współczynniki f; traktuje się jako niewiadome. Wstawiajac skończone sze- 


regi (28), (30) i (31) do równania (22) i uwzględniając ortonormalność 
układu funkcji w * ux, otrzymujemy: 


M 
= k & s 33 
9, 2 hs fa dla k-042,..., M. (33) 
Układ ten można łatwo rozwiązać numerycznie. Najpierw jednak należy 
obliczyć współczynniki gy i huk, co stanowi zasadniczy punkt powyższej 


metody. Zgodnie z równaniami (25) i (31) funkcja rozkładu energii %*(e) 
ma postać: 


* = z > ^ Y 
6*0 
gdzie f jest funkcja odwrotna do funkcji f, natomiast 


l * * 
kim 0c (5) =X (E) 
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Tak więc obliczenie funkcji y*(e) tą metodą sprowadza się do wykonania 
następujących operacji: 

1) obliczenia M + 1 calek (32), 

2) obliczenia (M +1)? calek (29), 

3) rozwiązania układu równań liniowych (33), 


4) obliczenia funkcji y*(e) danej wzorem (34). 

W pierwszych numerycznych obliczeniach wykonanych na bazie tej teorii 
autorzy przyjęli do opisu adsorpcji lokalnej równanie Langmuira (5) za- 
kłacając, że energia adsorpcji zmienia się w przedziale (0, oo). W nume- 
rycznych obliczeniach [38] wykorzystano dane eksperymentalne Draina 
i Morrisona [25], opisujące adsorpcję argonu na rutylu w temp. 85?K. 
Układ ten jest szczególnie odpowiedni do tego celu ze względu na dużą 
dokładność pomiarów i dużą liczbę punktów eksperymentalnych w szero- 
kim zakresie ciśnień. Drain i Morrison określili funkcję y*(e) dla 
tego układu, ekstrapolując wyniki otrzymane metodą przybliżenia kon- 
densacyjnego do temperatury zera kezwzględnego. Funkcję rozkładu 
energii y*(e) dla tego układu obliczyli również Dormant i Adamson 
[24], stosując swoją metodę sukcesywnych, graficznych aproksymacji. Ich 
wyniki zgadzają się dość dobrze z wynikami otrzymanymi wcześniej przez 
Draina i Morrisona [25]. Ponieważ do opisu adsorpcji lokalnej 
użyto równania Langmuira, dlatego też punkty eksperymentalne do obli- 
czeń brano z przedziału ciśnień poniżej 56 mm Hg. Przy tej wartości ciś- 
nienia ilość zaadsorbowana jest równa 755 cm?. Wartość ta zgadza się 
z pojemnością monowarstwy otrzymaną metodą BET przy przybliżonym 
założeniu, że powierzchnia rutylu jest homogeniczna. Funkcja rozkładu 
energii y*(e) znaleziona przez Rudzińskiego i współprace. [38] różni 
się nieco od funkcji rozkładu energii otrzymanej przez Dormanta 
i Adamsona [24] dla tego samego układu (patrz ryc. 3). Różnice te 
wynikają stąd, że do opisu adsorpcji lokalnej stosowano w obu przypad- 
kach różne izotermy adsorpcji. 

Następna seria obliczeń dotyczyła wpływu przedziału całkowania na 
obliczaną funkcję x*(e) [41]. W obliczeniach wykorzystano dane ekspery- 
mentalne Payne'a, Singa i Turka [42] dla układu adsorpcyjnego 
argon—żel krzemionkowy w temp. 77”K. Najniższą wartość energii ad- 
sorpcji oszacowano metodą przybliżenia kondensacyjnego, a stałą K zna- 
leziono z formuły Hobsona [43]. Z tak znalezionymi parametrami e, 
i K wykonano odpowiednie obliczenia numeryczne 4*(e) dla trzech różnych 
wartości eg. Wyniki tych obliczeń sugerują, że górna granica całkowania 
eg nie wpływa zasadniczo na kształt obliczonej funkcji x*(e). Z teoretycz- 
nych rozważań Heera [28] wynika, że na bardzo heterogenicznych po- 


184 W. Rudziński, J. Narkiewicz, M. Jaroniec 


20 25 ?B 32 36 
E kcal/mol 


Ryc. 3. Funkcje rozkładu energii adsorpcji dla argonu na rutylu w temp. 85°K; 
1 — funkcja z*(e) obliczona przy pomocy równania (34), 2 — funkcja Dormanta 
i Adamsona [24] 


wierzchniach, przy bardzo niskich ciśnieniach adsorbatu adsorpcja jest 
prawie skokowa. Przy umiarkowanie niskich p płaty z najwyższą energią 
adsorpcji są już zapełnione i dlatego nie wpływają na proces adsorpcji. 
Inaczej mówiąc, dane adsorpcyjne z obszaru średnich i niezbyt niskich 
ciśnień nie zależą od wartości maksymalnej energii adsorpcji eg. 

Z drugiej strony należałoby przypuszczać, że najniższa wartość energii 
adsorpcji e, będzie miała znaczny wpływ na pomiary adsorpcji przy śred- 
nich i niezbyt niskich ciśnieniach. Sokołowski i Jaroniec [41] 
potwierdzili to przypuszczenie wykonując modelowe obliczenia 4*(e) dla 
trzech różnych wartości e, przy ustalonym eg. Różnice pomiędzy funkcja- 
mi yx*(e), znalezionymi dla różnych wartości e, okazały się wyraźnie wi- 
doczne. 

Podsumowując, problem wpływu granic całkowania na rezultaty obli- 
czeń 4*(e) powinien być rozważany z uwzględnieniem wszystkich ,termo- 
dynamicznych warunków”, a więc i badanego przedziału ciśnień i tem- 
peratur. 

Przedstawione wyżej numeryczne metody obliczania wydają się bar- 
dzo użyteczne i stabilne, jeśli chodzi o błędy doświadczalne. Są one rów- 
nież zupełnie ogólne i matematycznie jasne, chociaż wymagają zaawanso- 
wanych i pochłaniających dużo czasu obliczeń numerycznych. 
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Ostatnio Jaroniec i Rudziński [44, 45] zaproponowali uprosz- 
czoną wersję tej metody, w której lokalną izotermę adsorpcji rozwija się 
w pojedynczy szereg Fouriera 


oo 


8, (p.F(6)) =) rle) "b [6] | (35) 


0*0 


Przypuśćmy, ze znamy wartości N;(p) dla ciągu wartości p;, j—0,1, ..., M. 
Podstawiając równania (31) i (35) do równania (19) otrzymamy układ rów- 
nań liniowych ze względu na nieznane współczynniki fi: 


M 


Nlp) =} rs(gi)y f, dla j=01...M (36) 
0:0 

W metodzie tej można użyć każdego ortonormalnego układu funkcji 
[ut] .Jaroniec i Rudziński [44] wykazali, ze w przypadku gdy 
©,(p, e) jest izotermą Langmuira lub BET najdogodniejszy jest zbiór funk- 
cji Lagrange'a [45]. W przedstawionej uproszczonej metodzie ilość ope- 
racji numerycznych znacznie się zmniejsza. Problem wyznaczania ;x*(e) 
sprowadza się do: 

1) obliczenia M--1 współczynników r; danych pojedynczymi całkami, 

2) rozwiązania układu równań liniowych (36), 

3) obliczenia funkcji x*(e) zgodnie z równaniem (34). 
Można łatwo zauważyć, że ta uproszczona metoda jest o wiele prostsza 
i mniej czasochłonna, a ponadto daje się matematycznie udowodnić, że 
przy M=oo obie te metody, zarówno oryginalna, jak i uproszczona, po- 
winny dawać dokładnie takie same rezultaty. Nawet jeśli M jest skończo- 
ne i niezbyt duże, wyniki otrzymane przy użyciu pełnej i przybliżonej 
wersji tej metody są dość zbliżone. 

Obszerne numeryczne badania Jarońca i Rudzińskiego [45] 
prowadzą do następujących wniosków: 


1. Pełmą wersję tej numerycznej metody należy stosować wtedy, gdy 
liczba punktów eksperymentalnych jest duża, ale pomiary są obciążone 
znacznymi błędami. 

2. Jeśli liczba punktów eksperymentalnych jest mała, ale nie spodzie- 
wamy się większych błędów eksperymentalnych, wtedy bez obawy może- 
my stosować uproszczoną wersję tej metody. 

Ilościowa charakterystyka efektów heterogeniczności jest jednym z 
ważnych problemów adsorpcji. Badania na ten temat są nadal prowadzo- 
ne. Najnowszymi proponowanymi metodami wyznaczania 4(:) są metody 
Rossa i Morrisona [46] oraz Landmana i Montroll [47]. 
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Pierwsza z nich jest w pewnym stopniu analogiczna do metody Van D o n- 
gena [27]. Druga zaś bazuje na procedurze rozwiązywania równania cał- 
kowego typu (1) proponowanej przez Wienera i Hopfa [48]. Oby- 
dwie metody nie znalazły jeszcze szerszego uznania i nie będą dyskuto- 
wane. 


PODSUMOWANIE 


Przedstawione wyżej numeryczne metody obliczania funkcji rozkładu 
energii x(e) z eksperymentalnycn danych adsorpcji stanowią poważny krok 
naprzód w ilościowych badaniach heterogeniczności powierzchni. Jak już 
wspomniano we wstępie, mogą one być wykorzystane do wyznaczania 
wielu innych wielkości fizykochemicznych, które są definiowane w spo- 
sób analogiczny do równania całkowego (1). Niektóre z dyskutowanych 
metod były szeroko stosowane w badaniach adsorpcji. Należą tutaj: me- 
toda Adamsona [18], Rossa [33], metoda Van Dongena zmo- 
dyfikowana przez Jarońca [30] oraz metody bazujące na rozwinięciach 
funkcji podcałkowych w szeregi Fouriera [38, 45]. Do najmniej praco- 
chłonnych metod należy metoda Rossa [33] oraz modyfikowana przez 
Jarońca metoda Van Dongena [530]. Najbardziej interesująco za- 
powiadają się metody bazujące na rozwinięciach funkcji podcałkowych 
w szeregi Fouriera. Tym niemniej wymagają one jeszcze wielu badań na- 
tury numerycznej, zmierzających do opracowania uniwersalnych proce- 
dur oraz szerszego wytestowania ich na większej liczbie układów adsorp- 
cyjnych. 

Dalsze badania nad tym problemem powinny zmierzać do opracowa- 
nia standartowych metod wyznaczania funkcji rozkładu energii adsorpcji. 
Zalążkiem takiej metody, służącej do charakterystyki gaussowskich ad- 
sorbentów, jest metoda Rossa. Następnym ważnym problemem jest po- 
równanie tego typu rozkładów z innego typu badaniami, np.: badaniami 
kalorymetrycznymi, rentgenograficznymi, spektroskopowymi itp. 

Interesująco zapowiada się porównanie funkcji rozkładu energii ad- 
sorpcji z rozkładami porów. Pozwolą one wyodrębnić efekty heterogenicz- 
ności związane ze strukturą geometryczną adsorbentu i chemiczną naturą 
powierzchni. Na obecnym, etapie badań nad rolą heterogeniczności w ad- 
sorpcji zachodzi potrzeba prowadzenia wnikliwych doświadczeń, pozwala- 
jących na odróżnienie efektów heterogeniczności od efektów związanych 
z oddziaływaniami w fazie adsorpcyjnej. Obydwa te efekty mogą się na- 
wzajem nakładać lub identyfikować. Tego typu badania są obecnie pro- 
wadzone w naszym Zakładzie. Jest rzeczą oczywistą, że właściwa ocena 
heterogeniczności powierzchni pozwoli na dokładniejszy opis szeregu zja- 
wisk zachodzących na granicy faz. 


wp 
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PE3IOME 


Obcy»XNAIOTCA BbIHMCJIMTeJbHble METOĄbI OripegeJieuMs (QYHKIMM pacrpe7eJeuus 


3Hepruu a4CcOopónun M3 3KCIepMMeHTAJIbHBIX U30TepM A4COPÓNNN. 


SUMMARY 


The critical review of numerical methods for evaluating the distribution functions 


of adsorption energy from experimental adsorption isotherms is given. 


